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1.  Introduction 
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In  diving  research  there  has  been  an  Interest  In  the  pulmonary  effects 
of  100%  0^  at  depth.  In  addition  to  pulmonary  effects  convulsions  also 
occur.  The  time  to  onset  of  convulsions  Is  reduced  as  depth  Is  Increased. 
In  this  report  we  develop  a mathematical  description  of  the  time  to  onset 
of  convulsions  based  on  the  data  of  Bradley  shown  In  Table  1.  Note  that 
the  time  to  onset  of  convullsons  is  given  for  each  of  23  dogs  at  a depth  of 
45  feet,  15  dogs  at  a depth  of  60  feet  and  10  dogs  at  a depth  of  86  feet. 

At  lower  depths  the  pulmonary  changes  develop  before  the  onset  of 
convulsions  and  may  actually  prevent  the  onset  of  convulsions.  This 
feature  has  not  yet  been  incorporated  into  the  model  except  that  we  do 
introduce  a form  for  the  depth  relationship  which  extrapolates  to  no 
convulsions  at  a zero  depth.  Several  generalizations  of  the  mathematical 
model  which  are  Introduced  provide  a basis  for  evaluations  of  the  model. 

2.  Descriptions  of  Time  to  Onset  of  Convulsions 

For  each  depth  we  assume  that  the  time  t to  onset  is  a random  variable 
which  has  a density  function  f(t).  The  density  function  f(t)  is  a 
mathematical  function  such  that  the  area  given  by 

f(x)dx 

ti 

is  the  probability  of  convulsions  In  the  interval  tj  to  t2.  We  may 


designate  this  probability  by  Pr(ti<t<t2). 
cumulative  distribution  function 

F(t)  = /*■  f(x)dx 

then 

Pr{ti<t<t2}  = F(t2)  “ F(ti) 
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and  F(t)  is  the  probability  of  onset  of  convulsions  before  t.  If  one  groups 
the  data  for  time  to  onset  of  convulsions  into  Intervals  the  relative 


frequency,  l.e.  the  number  of  animals  convulsed  between  tj  and  t2  divided 
by  the  number  of  animals  studied,  corresponds  to  the  observed  probability 
of  convulsions  in  the  interval.  Likewise  the  proportion  of  the  total 
number  of  animals  that  convulse  before  t corresponds  to  F(t). 

An  equivalent  alternative  description  Is  provided  by  the  risk  or 
hazard  function  for  convulsion 
h(t)  = f(t)/(l-F(t)) 

which  is  the  likelihood  of  the  onset  of  convulsions  at  t given  that 
convulsions  have  not  yet  occurred. 

Our  preliminary  perception  of  the  problem  was  that  the  risk  of 
convulsions  initially  would  be  at  or  near  zero  and  increase  with  time. 

Also  it  was  felt  that  onset  of  convulsions  could  be  thought  of  as  occurring 
as  the  result  of  a chain  one  or  more  reactions  or  processes.  From  this 
idea  we  formulate  a mathematical  form  which  describes  the  random  time  to 
onset  of  convulsions. 

3.  The  N-Reaction  Model 

Following  the  development  of  Bailey  et  al.  [1974]  we  equate  the  time 
to  onset  of  convulsions  with  the  transit  time  through  a chain  of  first 
order  consecutive  reactions  in  which  the  reaction  constants  are  taken 
to  be  equal.  The  resulting  cumulative  distribution  function  is 


n-1 

F (t)  - 1 - exp(-kt)  I — Tf — 

” i=0 

where  k is  the  reaction  constant  and  n is  the  number  of  reactions.  From 
this  we  find  the  density  function 
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d F_(t) 


consequently 


k (kt)“ 
(n-1) ! 


exp (-kt) 


. k (kt)"‘^  (kt)^  -1 

1=0 


For  n = 1,  hi(t)  = k.  That  Is  the  hazard  function  is  a constant.  For 

n = 2 or  greater  this  hazard  function  is  zero  at  t = 0 and  Increases  to 

the  asymptote  k as  t -»■  “>.  In  particular  for  n = 2 

2 

" 1 + kt 

and  for  n = 3 


h3(t)  - 


1 + kt  + (kt)2/2 


The  distribution  of  t is  further  characterized  by  either  the  moments  about 


the  origin 


, (v  + n-1) ! 
(n-1) ! k'' 


V = 1,  2,  3, . . . 


or  the  cumulants 


n(v-l) ! 


which  are  readily  obtained  from  the  moment  generating  function 


(k-s)“ 

We  use  these  expressions  to  obtain  the  mean  time  to  onset  of  convulsions 

> n 
Pi  - k 

and  the  variance  for  the  time  to  onset  of  convulsions 


Var(t)  - -7 
k 
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. Evaluation  of  the  N-Reaction  Model  for  the  Analysis  of  Time  to 


j 

* 

1 
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Onset  of  Convulsions 

Several  generalizations  of  the  n-reaction  model  were  needed  to 
evaluate  the  model.  First  we  added  a hazard  component  a to  the  risk 
function  for  n = 2 and  3.  So  that  for  t = 0 the  risk  would  be  a.  The 
estimate  of  a was  computationally  indistinguishable  from  zero  and  we 
concluded  that  the  constraint  that  the  model  have  an  initial  risk  of 
zero  in  these  cases  was  justified.  Furthermore  as  is  shown  in  Table  2 
the  best  fit  for  the  n-reactlon  model  is  for  n = 2. 

The  form 

k2  t 

"2(t)  = 


was  generalized  to  the  form 


V ^ (t/«) 

- 1 + (t/a)  • 

Furthermore,  the  Integral  of  h2{t)  designated  by 


H2(t)  = h2(t)dt 

= k [t  - a ln(l  + (t/a))] 

g 

was  generalized  by  replacing  (t/a)  by  (t/a)^  so 
H2(t)  = k[t^  - a^  In  (1  + (t/a)^)]. 

This  model  was  fitted  assuming  a common  S for  each  of  the  three  depths. 
The  log  likelihood  for  this  model  was  -285.98  which  is  only  a slight 
improvement  over  the  likelihoods  for  the  corresponding  (n=2)  2 and  3 
parameter  models  as  displayed  in  Table  2. 

5.  The  Depth  Relationship 

In  our  n-reaction  model  we  make  k a logistic  function  of  depth. 
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7.  Additional  Observations 
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Following  the  analysis  of  the  observations  in  Table  1,  an  additional  j 

seven  dogs  were  observed  at  a depth  of  30  feet.  There  were  two  dogs  * I 

observed  for  10  hours,  four  for  15  hours  and  one  for  18  hours.  In  each 
case  there  was  no  onset  of  convulsions.  Although  these  data  are  not  out 
of  line  with  the  predictions  given  in  Table  3,  estimates  of  the  para- 
meters for  our  2-reaction  model  were  reevaluated  using  the  additional 
data.  The  revised  parameter  estimates  are  0f  = 5.47  and  02  = 1.129  with 
the  approximate  covariance  matrix 
0f  02 

0.0384  -0.0254 

0.0175. 

Revised  predictions  are  shown  in  Table  4.  The  data  for  the  30  foot 
depth  is  in  a range  of  experimentation  which  is  more  difficult  because 
of  the  longer  time  to  onset  and  the  predicted  Increase  in  variability. 

Futhermore,  the  very  range  of  difficult  experimentation  is  one  in  which 
there  is  the  greatest  practical  interest. 
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Table  1.  Observed  Time  to  Convulsions  for 


n 

Average 


each  dog  at  each  of  three  depths 
Time  to  convulsion  (hours,  minutes) 


45  feet 

60  feet 

86  feet 

0:28 

0:51 

0:48 

1:05 

5:31 

0:27 

6:38 

5:02 

0:25 

0:26 

0:49 

1:52 

5:15 

2:40 

0:24 

8:05 

2:15 

1:06 

8:39 

2:01 

1:28 

1:17 

4:42 

0:17 

1:00 

5:20 

1:22 

6:30 

3:04 

0:46 

1:00 

4:25 

4:38 

1:32 

1:00 

0:49 

7:46 

0:54 

7:19 

0:26 

5:38 

4:51 

7:26 

7:02 

3:28 

1:10 

4:13 

5:21 

23  15  10 

4:22  2:41  0:54 


Table  2.  Log  likelihood  for  the  n-reaction  model  in  which 
it  shown  that  n = 2 gives  the  best  fit  and  that  depths 


( 


are  significantly  different. 


n 

depths  pooled 

1 parameter 

depth  fitted  to 
the  logistic  form 
2 parameters 

Log  likelihood 
depths  treated 
separtely 

3 parameters 

1 

-294.77 

>-292.09* 

-292.09 

2 

-293.12 

-288.24 

-287.71 

3 

-299.67 

-292.30 

-291.49 

*The 

form  of  the  2 parameter 

model  was  not  run 

since  it  can  not  give  a 

lower  value  than  obtained  for  the  3 parameter  model  and  the  maximized 
likelihood  is  clearly  at  n = 2. 
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Table  3.  Predictions  at  selected  depths  for  the  upper  asymptote  k 
of  the  risk  function,  the  average  time  to  convulsions  and  the 
variance  of  the  time  to  convulsions  based  on  our  maximum 


likelihood  estimates  for 

the  2-reaction 

model. 

Depth,  in 

feet  k(min 

U (min) 

Var (min^) 

5 

3.79x10"^ 

5.28x10^ 

1.39x10^ 

10 

1.99x10"^ 

1.00x10^ 

5.05x10^ 

15 

5.25x10"^ 

3.81x10^ 

7.25x10^ 

20 

1.05x10"^ 

1.91x10^ 

1.83x10^ 

30 

2.76x10"^ 

7.26x10^ 

2.63x10^ 

45* 

7.24x10"^ 

2.76x10^ 

3.81x^.^ 

60* 

1.43x10"^ 

1.40x10^ 

9.75x10^ 

86* 

3.33x10"^ 

6.00x10 

1.81x10^ 

100 

4.70x10"^ 

4.30x10 

9.04x10^ 

200 

2.06x10“^ 

9.70 

4.71x10 

500 

7.00x10"^ 

2.86 

4.00 

1000 

9.24x10"^ 

2.16 

2.34 

*Indlcates 

observed  depths  as  shown  in 

Table  1. 
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Table  4.  Revised  predictions  based  on  the  additional  seven 
observations  at  30  feet.  For  selected  depths  the  estimate 
of  the  upper  asymptote  k of  the  risk  function,  the  average 
time  to  convulsions  and  the  variance  of  the  time  to 
convulsions  based  on  our  revised  maximum  likelihood 
estimates  for  the  2-reaction  model. 


Depth,  in  feet 

k(mln 

u(mln) 

Var (min^) 

5 

6.47x10”^ 

3.09x10^ 

4.78x10^° 

10 

5.52x10*^ 

3.62x10^ 

6.56x10® 

15 

1.94x10"^ 

1.03x10^ 

5.34x10^ 

20 

4.71x10"^ 

4.24x10^ 

9.01x10® 

30 

1.65x10"^ 

1.21x10^ 

7.35x10® 

45 

5.76x10"^ 

3.47x10^ 

6.03x10® 

60 

1.39x10"^ 

1.44x10^ 

1.03x10^ 

86 

4.12x10“^ 

4.85x10 

1.18x10® 

100 

6.41x10"^ 

3.12x10 

4.86x10^ 

200 

3.69x10“^ 

5.42 

1.47x10 

500 

9.09x10“^ 

2.20 

2.42 

1000 

9.88x10"^ 

2.02 

2.05 

! 


onset  of  convulsions 
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